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ABSTRACT 

We have discovered a triple-peaked X-ray burst from the low-mass X-ray binary 
(LMXB) 4U 1636-53 with the Rossi X-ray Timing Explorer (RXTE). This is the first 
triple-peaked burst reported from any LMXB using RXTE, and it is only the second 
burst of this kind observed from any source. (The previous one was also from 4U 
1636-53, and was observed with EXOSAT.) From fits to time-resolved spectra, we 
find that this is not a radius-expansion burst, and that the same triple-peaked pattern 
seen in the X-ray light curve is also present in the bolometric light curve of the burst. 
Similar to what was previously observed in double-peaked bursts from this source, 
the radius of the emitting area increases steadily during the burst, with short periods 
in between during which the radius remains more or less constant. The temperature 
first increases steeply, and then decreases across the burst also showing three peaks. 
The first and last peak in the temperature profile occur, respectively, significantly 
before and significantly after the first and last peaks in the X-ray and bolometric light 
curves. We found no significant oscillations during this burst. This triple-peaked burst, 
as well as the one observed with EXOSAT and the double-peak bursts in this source, 
all took place when 4U 1636-53 occupied a relatively narrow region in the colour- 
colour diagram, corresponding to a relatively high (inferred) mass-accretion rate. No 
model presently available is able to explain the multiple-peaked bursts. 

Key words: stars: neutron — X-rays: binaries — X-rays: bursts — stars: individual: 
4U 1636-53 



1 INTRODUCTION 

Thermonuclear, type-I, X-ray bursts (also known as 
X-ray bursts, or simply bursts), are due to unsta- 
ble burning of H and He on the surface of accreting 
neutron stars in low-mass X-ray binaries (LMXBs). 
During these bursts, the X-ray flux increases by 
factors of ^^10-100 within a second or less, after 
which it returns back to its pre-burst level over a 
time-span of tens of seconds, fol lowing an aj^prox- 
imately exponential decay (e.g., Lewin et al. 1993|; 
[Strohmayer &: Bildsten 2003; Galloway et al. 20oi). 

One of the best studied sources of X-ray bursts is the 
LMXB 4U 1636-53. Also kno wn as V801 Ara, 4 U 1636-53 
was discovered with OSO-8 l| Swank et al.l [l976h . and was 



subsequently studied in great detail using obser vations with 
SAS- 3, Hakucho, Tenma, and EXOSAT (see iLewin et al] 
Il987l . for a review). The orbital per iod of this binary system 
is 3.8 hr (Ivan Paradiis et al] ll990'). and the spin period of 
the neutron st ar is 581 Hz |S trohmavcr et al. 1998a. b). Us- 
ing EXOSAT, iDamen et 311(119891 ') detected 60 bursts from 
this source between 1983 and 1986; from observations with 
the Rossi X-ray Timing Ex plorer (RXTE), 172 b ursts were 
detected up to 2007 June 3 (|Gallowav et al.l l200^. and more 
than 250 bursts including data taken after that date (Zhang 
et al., in prep.). Most of these X-ray bursts have standard, 
single-peaked, fast rising and exponentially decaying light 
curves. Three superbursts have also been detected i n 
this source d Wiinands 200lJ; iKuulkers et al.l [2004'). 
iLinares et al.i (2009) detected a 5.4 minutes recur- 
rence time between two X-ray bursts in 4U 1636— 
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53, which is the shortest observed recurrence time 
in any neutron star low-mass X- ray binary. 

However, there have also been several detections of 
bursts from 4U 1636—53 that have double-p eaked light 
curves, both with EXQSAT (ISztaino et al.lll985l) and RXTE 
ijBhattacharvva fc StrohmaveJ l2006bl : iGallowav et al.l 
l2008t ). Only one triple-peaked X-ray burst has ev er been 
detected from this source by Ivan Paradiis et al] (|l986l ) 
using EXOSAT observations. Several models have been 
proposed to explain the doubl e -peake d bursts. For instance, 
iBhattacharvva fc Strohmaverl (|2006bl ) suggested that these 
bursts are due to stalling of the nuclear-burning front near 
the equator of the neutron star following igniti on near the 
neutron-star pole, whereas iFisker et al.l (|2004h suggested 
that these bursts are caused by a waiting point in the 
thermonuclear reaction chain. 

In this paper, we report the first detection of a triple- 
peaked X-ray burst in 4U 1636—53 with RXTE, and only 
the second ever detected from this source. The paper is or- 
ganized as follows: We describe the observations and data 
analysis in ^ we show the results in i|3] and finally we dis- 
cuss our findings in ^ 



2 OBSERVATIONS AND DATA ANALYSIS 

We analysed all data available from the RXTE Proportional 
Counter Array (PCA) of 4U 1636-53 as of March 15 2009. 
The PCA consists of an array of five coUimated proportional 
counter units (PCUs) operating in the 2—60 keV range, with 
a total effective area of a pproximately 6500 c m^ and a field 
of view of ~ 1° FWHM jjahoda et al.lll996l ). For each ob- 
servation we calculated X-ray colours from the Standard2 
data (16-s time-resolution and 129 energy channels in the 
2 — 60 keV band), and produced 1-s light curves from the 
Standardl data (0.125-s time resolution with no energy res- 
olution). To calculate the colours of the source, we first ex- 
amined the 1-s Standardl light curves to identify and re- 
move any X-ray burst from the data, and we subsequently 
computed light curves in four energy bands every 16 s from 
the Standard2 data, separately for each PCU detector that 
was active during an observation. We subtracted the back- 
ground contribution from each light curve, and we corrected 
the count rates for dead time. We defined the soft colour as 
the count rate in the 3.5—6.0 keV band divided by the count 
rate in the 2.0—3.5 keV band, and the hard colour as the 
count rate in the 9.7—16.0 keV band divided by the count 
rate in the 6.0—9.7 keV band. We used data of the Crab pul- 
sar and nebula taken close to each observ ation of 4U 1636- 
53 to correct for instrumental effects (e.g. lAltamirano et al.l 
|2008| ). All colours of 4U 1636-53 presented in this paper are 
therefore normalised to the colours of Crab. Whenever we 
found an X-ray burst, we used 80 s of persistent emission 
just before the burst to represent the colours of the source 
at the time the burst started. 

In this paper we concentrate on an unusual X-ray burst 
that took place at 18:23:25 UTC on December 11 2006 (MJD 
54080.76626; ObsID 92023-01-44-10). To calculate the bolo- 
metric flux of the persistent emission before the burst, Fpboi, 
we extracted energy spectra of ~ 800 s before the onset of 
the burst from the two main instruments on board RXTE, 
the PCA and the High Energy X-ray Timing Experiment 



(HEXTE). HEXTE consists of two clusters of four Nal/Csl 
phoswich scintillation detec tors that are sensitiv e to X-rays 
in the 15 - 250 keV range (|Gruber et al.lll996l ). Cluster A 
of HEXTE stopped working properly in October 2006, an 
therefore only cluster B was available for our analysis. For 
the PCA spectrum we analysed the Standard2 data of PCU 
2, while for the HEXTE spectrum we analysed the Standard 
data of cluster-B. We extracted PCA and HEXTE back- 
grounds, and produced response matrices for both instru- 
ments following the instructions in the RXTE pages. We 
added a 0.6% systematic error to the PCU-2 spectrum, but 
no systematic error was applied to the HEXTE spectrum. 

To study this burst in detail, we analysed the PCA 
Event data, E_125us_64M_0_ls, in which each individual 
photon is time tagged at a ~ 122 fis time resolution in 64 
energy channels between 2 — 60 keV. We used all PCUs 
that were operating at the time of the X-ray burst (PCUs 
0,1,2 and 4) to produce 0.25-s resolution light curves in the 
full PCA band and in two other bands, 2.0 — 3.5 keV and 
6.0 — 9.7 keV. For the time-resolved spectral analysis of this 
burst we extracted spectra in 64 channels every 0.25 s from 
the Event data of all available PCUs. Since the light curve 
of the decay of the burst is quite smooth, we extracted spec- 
tra over somewhat longer intervals in the tail of the burst 
to compensate for the lower count rates. We used a time 
resolution of 0.5 s, 1 s and 2 s, respectively, for each of the 
following time intervals: 20.5 — 25.5 s, 25.5 — 28.5 s, and 
28.5 — 32.5 s after the start of the burst. We generated an 
instrument response matrix, and we fltted the spectra us- 
ing XSPEC version 12.4.0. Because of the short exposure, 
in this case the statistical errors dominate, and therefore 
we did not add any systematic error to the spectra. We re- 
stricted the spectral fits to the energy range 3.0 — 20.0 keV. 
During the ~ 800 s previous to the onset of the burst, the 
light curve and colours of 4U 1636-53 were consistent with 
being constant, therefore we extracted a spectrum of those 
data to use as background in our fits; this approach is well 
establis hed as a standard pr ocedure in X-ray burst analy- 
sis (e.g. iKuulkers et al.ll2003 ). We note that this procedure 
fails if the blackbody emission during the burst comes from 
the same source that produces the blackbody emission seen 
in the persistent emission, since the difference between two 
black body spectra is not a blackbody l|van Paradiis fc LewinI 
Il986l ). This effect is significant only when the net burst emis- 
sion is small, and therefore the problems may arise only at 
the start and tail of the burst, when the burst emission is 
comparable to the p ersistent emission (see the discussion in 
iKuulkers et al.ll2002l ). We corrected each spectrum for dead 
time using the methods supplied by the RXTE team. 

We computed Fourier power density spectra for the du- 
ration of the burst using the event data over the full PCA 
band pass. We calculated the power spectra from 1-s data 
segments, setting the start time of each segment to 0.125 
s after the start time of the previous segment. Note that 
because of this, the individual power spectra are not inde- 
pendent. 

3 RESULTS 

In the top panel of Figure [T] we show the light curve of 
the burst at a resolution of 0.25 s. The light curve dis- 
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Figure 1. Top panel: Light curve in tiio 2—60 keV range of the 
triple-peaked burst observed in 4U 1636—53. Time is given in 
seconds from 2006 December 11 at 18:23:25 UTC. Bottom panel: 
The time evolution of the hardness ratio during the triple-peaked 
burst in 4U 1636—53. The solid vertical lines indicate the times of 
local maxima in the light curve, extending into the hardness-curve 
panel. The dashed vertical lines indicate the times of local minima 
in the light curve, extending into the hardness-curve panel. The 
left-most dotted-dashed vertical line indicates the time of the first 
maximum in the hardness curve extended into the light-curve 
panel. We also plotted the Itr error bars. 



plays three peaks, similar to wh at was observed dur i ng the 
triple-peaked burst reported by Ivan Paradiis et all (|l986l ) 
from EXOSAT data of this source, and we therefore call 
this a triple-peaked burst. Out of the three peaks in the 
burst, the first peak is the brightest in the full PGA band 
(~ 2000 counts s~^ PCU~^) and the second is the weak- 
est (~ 1700 counts s~^ PCU~^). The separation in time 
between the first and third peak in the triple-peaked burst 
re ported here is ^ 8 Sj_whereas in the triple-peaked burst 
in Ivan Paradiis et al.l (| 19861" ) it twice as long, ~ 17 s. In the 
bottom panel of Figure [T] we show the hardness curve during 
the burst. The hardness is defined as the count-rate ratio in 
the 6.0 - 9.7 keV and 2.0 - 3.5 keV energy bands. 

Besides the main three peaks, both curves show struc- 
ture on shorter time scales, which makes it difficult to ac- 
curately identify the times of the maxima of each peak. It 
is nevertheless apparent from this Figure that the first peak 
in the light curve occurs significantly later than the first 
peak in the hardness curve. Since the hardness is a mea- 
sure of the temperature of the emitting surface (see below) , 
this means that the emitting surface reaches the maximum 
temperature before the burst is the brightest. The first peak 
in the light curve occurs about 2.5 s after the first peak in 
the hardness curve. A similar delay between light curve and 
hardness curve h as been observed in phot ospheric radius- 
expansion bursts jStrohmaver et aL I ll998bl ), and is indica- 
tive of the start of the radius-expansion phase but, as we 
discuss below, this triple-peaked burst is not a radius ex- 
pansion burst. After this initial mismatch between light and 
hardness curves, the other two maxima and the two minima 
of the light curve appear to coincide with the corresponding 
maxima and minima in the hardness curve. 



3.1 Pre-burst spectrum 

We fitted the PGA (in the 3.0-25.0 keV range) and HEXTE- 
B (in the 15 — 200 keV range) pre-burst spectra simul- 
taneously with a model consisting of a combination of a 
blackbody, a power law, and a Gaussian emission line fixed 
at 6.4 keV, all affected by interstellar absorption. We in- 
cluded a multiplicative factor in the model to account for 
possible systematic differences in the relative calibration 
of the effective areas of the two instruments. We found a 
good fit to the spectra, = 92 for 89 degrees of freedom 
(dof), using a model consisting of a blackbody with tem- 
perature kT — 1.92 ± 0.05 keV and bolometric fiux F^h = 
6.3 ± 0.7 X 10^^" erg cm^^ s^^, a power law with photon 
index F = 2.96 ± 0.05 and normafization TVpi = 1.77 ± 0.14 
photons cm^^ s^^ keV^^ at 1 keV, and a Gaussian line 
with a = 1.08 ± 0.18 keV and normalization A^g = 4.2 ± 
1.7 X 10~^ photon cm~^ s~^. We also included the effect 
of interstellar absorption using the cross-sections of 
iBalucinska-Churc h &: McCammon (1992) and solar 
abundances from ,Anders &: Gr evesse ( 1989, ) , with a 
hydrogen e quivalent column d ensity fixed at A'h = 0.36 x 
lO^'^ cm~^ (jPandel et al.ll2008). Th ese parameters are com- 
parable to those reported by Gackett et aL (20 0^ f or thi s 
source in the soft state. Following iGallowav et al.l l|2008t ). 
we used the unabsorbed 2.5 — 25 keV fiux, ^2.5-25 = 
1.68 ± 0.02 x 10^^ erg cm"'^ s~^, and a bolometric correc- 
tion Cboi ~ 1.38, to estimate the bolometric flux before the 
burst, Fpboi = 2.32 ± 0.03 x 10"^ erg cm"^ s"^ We 
note that the error o f Cboi can be as large as 40 % 
(| Galloway et al . 1 120081 ') ■ Using the Eddington lumi- 
nosity for a 1.4-Mq neutron star with a 10-km ra- 
dius, 7 = -Fpboi/-FEdd = 3.7 ±0.2 X 10~^ - 6.2 ±0.2 x 10"^ 
for a range of hydrogen mass fraction X = — 0.7. 



3.2 Time-resolved spectra during the burst 

We fitted the time-resolved net burst spectra with a single- 
temperature blackbody model (bbodyrad in XSPEG), as 
generally burst spectra are well fit by a blackbody. During 
the fitting, we kept t he hydrogen colum n density A^h fixed 
at 0.36 X lO^^cm"^ l|Pandel et al.l 120081 ). and to calculate 
the radius of the emitting area in k m, we assumed a dis- 
tance of 5.9 kpc l|Fiocchi et a l."2006Y The model provides 
the blackbody colour temperature (Tc) and a normalization 
proportional to the square of the blackbody radius (-Rbb) 
of the burst emission surface, and allows us to estimate the 
bolometric fiux as a f unction of time . 

It is well kno wn (|London et aLlll986l : [Titarchuklll994l : 
iMadei et aLll2004h that the colour temperature Tc obtained 
from fitting the continuum spectra of an X-ray burst is 
higher than the effective temperature Tctt that enters in the 
calculation of the bolometric flux. In order to get the ef- 
fective temperature Toff, we assumed a neutron star with 
a mass of I.4M0 and a radius of 10 km. Interpolatin g the 
ratio of Tc/T^a from the table of lMadei etai] l|2004l ). and 
using the bolometric flux from the fit, we calculated the "cor- 
rected" blackbody radius (see also van Paradiis ct aL , 19861 ). 
We found correction factors between 1.3 and 1.5. 

In Figure [5] we show the time evolution of the 
best-fitting spectral parameters: Effective temperature, 
Tcft , "corrected" blackbody radius, Rbh, and bolometric flux. 
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Figure 2. Time evolution of spectral parameters during the 
triple-peaked burst in 4U 1636-53. From top to bottom the panels 
show the effective temperature, the 'corrected' (see text for de- 
tails) radius of the emitting area, the bolometric flux, the 2 — 60 
keV light curve, and the reduced of the fits (for 18 dof). The 
time resolution of this plot is 0.25 s for times between s and 
20.5 s, 0.5 s for times between 20.5 s and 25.5 s, 1 s for times 
between 25.5 s and 28.5 s, and 2 s from 28.5 s until the end. The 
solid vertical lines indicate the times of local maxima in the X-ray 
light curve (second panel from the bottom), extending into the 
other panels. The dashed vertical lines indicate the times of local 
minima in the X-ray light curve extending into the other panels. 
The left-most and right-most dotted-dashed vertical line indicate 
the times of the first and last maximum in the temperature curve 
(upper panel) extended into the other panels. We also plotted in 
each panel the Icr error bars. 



Fboi- For reference we also plot the 2 — 60 keV light curve 
in this Figure. X-ray hght curves with multiple peaks are 
often a signature of radius expansion bursts, since during 
the radius expansion phase the temperature of the emit- 
ting surface decreases, and part of the rad iation falls outside 
the X-ray band (e.g., iLewin et al.lll993l ). In those bursts, 
however, the bolometric light curve shows a single peak. 
In this burst, the bolometric light curve also shows three 
peaks, at the same positions as the peaks in the X-ray 
light curve. This shows that the multiple peaks in the X- 
ray light curve are not due to photospheric radius expan- 
sion during the burst, which is consistent with the fact that 
the flux at the peak of this burst is about a factor 2 — 3 
smaller than in re al radius expansion bursts in this source 
l|Maurer fc Watts! 2008). The total burst fluence was 
Eh = 2.9 ± 0.1 X 10"^ erg cm^ . From this, and the peak 
of the bolometric flux, the characteristic time scale of 
the burst was r — -Bb/fpoak ~ 13.8 ± 0.7 s. Most of the 
fits are good, although there are a few cases with relatively 



large reduced values. (Only in six occasions the reduced 
is larger than 2.) We note that we did not add a sys- 
tematic error to the data (see SJJ]). We confirmed that the 
distribution of values is consistent with a x^ distribution 
with 18 dof, which is the number of degrees of freedom in 
the individual fits. 

The top panel of Figure [2] shows that, after an initial 
rise, on average the effective temperature decreases with 
time. The temperature profile shows three maxima, and as 
in the case of the hardness curve, the first peak of the ef- 
fective temperature takes place before the first peak of the 
X-ray and bolometric light curves. The second peak of the 
effective temperature profile more or less coincides with the 
second peak in the X-ray and bolometric light curves, while 
the two minima in the temperature profile also coincide with 
the times of minima in the bolometric and the X-ray light 
curves. The last peak in the temperature curve occurs sig- 
nificantly after the last peak in the light curve. 

The second panel of Figure [2] shows the evolution of 
the radius of the emitting surface, corrected for the fact 
that the spectrum during the burst is non-Planckian (see 
above) . On average, the radius increases steadily during the 
burs t, similar to what was observed in dou ble-peaked bursts 
(e.g. iBhattacharvva fc Strohmaveill2006al lbl). After initially 
rising, the radius stops increasing at the time of the first 
peak in the light curve, although there is no apparent change 
in the trend of the temperature at that point in time. For 
2 — 3 s the radius remains more or less constant, and after the 
second maximum in the light curve, it increases again. More 
or less at the time of the second minimum in the light curve 
the radius stops increasing, or slightly decreases, and after 
the third maximum the radius continues increasing until the 
end of the burst. (We note that the slight increase in the 
radius at the time of minimum temperature could be due to 
limitations in the colour correction factor.) 

3.3 Timing analysis of the burst 

We searched this burst for oscillations around the known 
spin frequency of the ne utron star in 4U 1636-53 (581 Hz; 
IStrohmaver et ahl Il998al ) . For this, we analysed the data 
within 15 s of the onset of the burst, and we searched only 
the frequency range 579 — 583 Hz. We chose a detection limit 
given by a single trial probability of 3 x 10~^ or less of hav- 
ing a given power due to random fluctuations in the data. 
This is equivalent to a ~ 4a detection limit considering the 
number of trials involved. We did not flnd any significant 
oscillations in this frequency range during the burst. The 
95% confidence upper limits to any oscillation in the fre- 
quency range 579 Hz to 583 Hz were 7.5%, 12.3% and 7.5% 
rms for each of the three peaks of the burst, respectively, 
and 12.2% and 11.3% rms for the rise and the decay of the 
burst, respectively. Previously, burst oscillations with am- 
plitudes between ^ 2% and ~ 9% were detected from 4U 
1636-53 l|Zhang et al.lll997l ). 

Most burst oscillations in 4U 1636-53 were so far 
dete cted in br ight, photospheric radius expansion bursts 
(Mau rer fc Wa tts 2008), with bolometric peak fluxes larger 
than ~ 5.0 x 10^** erg s~^ cm~^. One of the few excep- 
tions to this trend is the possible detection of oscillations 
during the firs t peak of a non-radius expansion d ouble- 
peaked burst in IBhattacharvva fc Strohmaveij (|2006lj ): The 
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bolometric flux at the peak in which osciUations were re- 
ported was ~ 1.6 X 10~* erg s~^ cm~^. The bolometric flux 
at the peak of the triple-peaked burst reported here was 
~ 2.1 X 10~* e rg cm~^, similar to that o f the do uble- 
peaked burst of iBhattacharvva fc Strohmaverl (|2006bl ') . and 
a factor 2 — 3 weaker than other photospheric radius ex- 
pansion bursts with oscillations in this source. The upper 
limit to the rms amplitude we found during the first maxi- 
mum of the triple-peaked burst of 4U 1636-53 reported here, 
7.5%, is lower than the rms amplitude of the oscillations 
during the first maximum of the double-peaked bu rst, 8.2%, 
reported bv lBhattacharvva fc Strohmaveij l|2006bl ) from this 
same source. 



3.4 Colour-colour diagram 

In Figure |3] we show the colour-colour diagram of all the 
RXTE observations of 4U 1636-53 to date (grey points), 
excluding times of X-ray bursts; the black crosses in the 
plot are the colours of the source at the start of an X- 
ray burst. There are more than 250 X-ray bursts shown 
in this Figure, from RXTE data taken before March 2009. 
The colour-colour diagram of 4U 1636-53 is typical of 
that of a low-luminosity LMX B, a so-called Atoll source 
iHasinger fc van der KM1989I ). The ~ 250 X-ray bursts ob- 
served with RXTE distribute more or les s uniforml y across 
the c olour-colour diagram (Belloni et al. 112007; M uno et al.l 
|2004| ). The triple-peaked burst, which we indicated using 
a open circle, is located close to the vertex of this colour- 
colour diagram, in an area where there are several other 
bursts. For comparison, we plot the approximate position of 
the triple-peaked burst d etected with EXOSAT in Figure [3] 
l|van der Khs et al.llT990l) . 



3.5 Bursts properties 

In order to compare the properties of single-peaked 
bursts and muti-peaked bursts that occurred in the 
same area of colour-colour diagram, we selected 
all the bursts which have hard color less than 0.7 
(see Figure [3|), and calculated their fluences, peak 
fluxes and peak temperatures. We found that: (i) 
Single-peaked bursts have fluences in the range 
2 X 10"* to 70 X 10"* ergs cm~^, whereas all the 
multi-peaked bursts have fluences that span a nar- 
rower range, from 23 x 10"* to 38 x 10"* ergs cm"^; 
(ii) single-peaked bursts show peak color tempera- 
tures from 1.3 to 2.7 keV, while the multi-peaked 
bursts have temperatures in the range 1.7 to 2.1 
keV, with the exception of the double-peaked burst 
from IBhattacharvva &: Strohmaverl (|2006bl ). which 
has a peak temperature of 2.37±0.04 keV. (iii) single- 
peaked bursts show peak fluxes between 3 x 10"^ 
and 82 x 10"^ ergs cm"^ s"^, while the multi-peaked 
bursts have peak fluxes from 21 x lO"'' and 28 x 10~' 
ergs cm~'^ s ~^, except for the same doubl e- peaked 
burst from IBhattacharvva 8z Strohmaverl (|2006bl ) 
which has a peak flux of 53 ± 4 x 10~® ergs cm~'^ s~^ 
(Zhang et al. in prep.). 



4 DISCUSSION 

We have detected a triple-peaked X-ray burst in the LMXB 
4U 1636-53. This is only the second triple-peaked burst ever 
detected from any LMXB, and the first one observed with 
RXTE. (The only other triple-peaked burst reported is also 
from 4U 1636-5 3, and was observed more than 20 years ago 
with EXOSAT: [van Paradiis et all [l983) . The bolometric 
light curve of the triple-peaked burst shows a profile similar 
to that of the X-ray light curve, which indicates that this 
was not a radius expansion burst (see ii3.2[) . This is consis- 
tent with the fact that the bolometric flux at the peak of 
this burst was a factor of 2 — 3 smaller than the bolometric 
flux observed from radius-ex pansion bursts from this source 
(e.g. iMaurer fc Wattsll2008l ). 

The temperature curve during this burst shows the 
same pattern of three maxima seen in the X-ray light curve, 
however the flrst and last peaks in the temperature curve oc- 
cur signiflcantly before and signiflcantly after the first and 
last peaks in the X-ray light curve, respectively. The radius 
of the emitting surface increases more or less steadily during 
the burst, with small intervals during which the radius re- 
mains approximately constant or perhaps decreases slightly. 

We found no oscillations in this burst, with upper lim- 
its that are consistent with oscillations detected in previous 
bursts in 4U 1636-53, but smaller than the amplitude of the 
oscillations see n during the non-radius expansion, do uble- 
peaked burst in IBhattacharvva fc Strohmaverl (|2006bl ). 

The triple-peaked X-ray burst took place when the 
source was at an intermediate position in the colour-colour 
diagram. There are several other X-ray bursts around this 
location in the colour-colour diagram, none of them showing 
a triple-peaked profile. If the position of the source along the 
C-like shape of the colour-colour diagram (Fig. |3]) is a mea- 
sure of mass accretion rate (with mass accretion rate increas- 
ing as the source moves along the branches from the upper- 
right to the l ower-left and then to the lower -right corner of 
this diagram; iHasinger fc van der Klislll989l l. this indicates 
that the triple-peaked nature of this burst is not associated 
to a rate of mass accretion onto the neutron star that is dif- 
ferent from that in other, single-peaked, bursts. The previ- 
ously detected triple-peaked burst in 4U 1636-53 took place 
when the source was located in the middle of the branch 
that extends to the right of th e position of the triple burst in 
this diagram (|van der Klis et a l. 1990) , when mass accretion 
rate was presumably higher. Also in the case of the EXOSAT 
observations, there are other n ormal bursts arouiid the po- 
sition of triple-peaked burst of Ivan Paradiis et al.l (|l986l ) in 
the colour-colour diagram. 

The triple-peaked burst started ~ 800 s from the begin- 
ning of the RXTE observation; from the absence of another 
X-ray burst in those ~ 800 s, the wait time from the previ- 
ous X-ray burst was > 800 s. The previous X-ray burst 
recorded by RXTE before the triple-peaked burst took place 
at 02:42:12 UTC on November 19 2006, therefore the wait 
time before the triple-peaked burst w as ^, 2 x 10^ s. Give n 
the typical burst rate in 4U 1636-53 (|GaUowav et al.ll2008l '). 
most likely several X-ray bursts were missed in between 
these two observations, and therefore we can only give a 
lower limit of 6.7±0.1 to the value of q = (-F'pboifw)/i5b 
of the triple-peaked burst. 

iBoirin et al] l|2007l ) reported X-ray bursts that occur in 
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Figure 3. Colour-colour diagram of all observations of 4U 1636—53 up to March 2009. The grey points represent the colours of the source 
from all RXTE observations available, excluding the times of X-ray bursts (see text). Each point in this diagram represents 256 s of data. 
We plot the representative Irr error bar in the upper left part of this Figure. The black crosses represent the colours of the persistent 
emission of the source at the onset of an X-ray burst. The open squares indicate double-peaked bursts (Galloway et al. 2008; Zhang et 
al. in prep.). The open circle indicates the triplc-pcakcd burst. The op en triangle shows the approximate location in this diagram of the 
triple-peaked burst observed with EXOSAT (,van dcr Klis et al.lll99oli . 



triplets; we note that theirs is probably a completely differ- 
ent phenomenon than the tri ple-peaked burst tha t we report 
here. The bursts detected bv lBoiriri et all l|200'i1 ) are sinele- 
peaked, with waiting times of ^ 12 min between the three 
components of the triplet. The burst we report here shows a 
triple-peaked structure within a single burst, with the time 
between peaks being less than ~ 10 seconds. 

Most type-I X-ray bursts display a single peak in 
the bolometric light curve, with a profile that shows a 
fast rise and an exponential decay (jCallowav et al.ll2008h . 
About a dozen double-peaked and two triple-peaked X-ray 
bursts have been detected from the LMXB 4U 1636—53 



llSztaino et al. 19851; 


Lewin et al.l 19871; van Paradiis et al. 


19861; iGallowav et al. 


I2OO8I; iBhattacharvva & Strohmaver 


2006allbl: Maurer & Wattd 20081, this oaDer). Double-Deaked 



bursts were observed fro m other sources besi des 4U 1636- 
53, e.g., 4U 1608-52 l|Penninx et al.1 1 19891 ) GX 17+2 
jKuulkers et al.1 |2002| ). and 4U 1709-267 (Ijonker et al.1 
I2OO4I ); but so far there is no report of a triple-peaked burst 
from another source besides 4U 1636-52. 



In the past 20 years, several models have been pro- 
posed to explain the doubl e-peaked non-radius expansion 
bursts. i Regev fc Livid (Il984l ) proposed a model in which the 
neutron-star surface is divided in two zones, such that the 
burst starts in one zone and moves into the other. In this 
model, a temperature gradient develops between the zone 
where the burst ignites and the zone onto which the burst 
flame expands. A heat transport impedance between the two 



zones causes a dip in the temperature and the light curve , 
and produces a double-peaked burst. iFisker et al.l (|2004l ) 
repeated these calculations using up two 200 zones, and 
found that the effect dis appears as the numb er of zones 
increases beyond ~ 25. iFuiimoto et al.l l|l988h explained 
the double-peaked bursts as a two-step energy generation 
due to shear inst abilities in the fuel on the stellar sur- 
face. Bhattacharvv a fc Strohmaverl (|2006al ) argued that this 
model cannot reproduce the double-peaked profiles that are 
observed, and that it is difficult to maintain sufficient un- 
burnt material in a layer above the burning flame, as re- 
quired for the instability to take place, without having the 
hot fuel mixing with the cold one. .Meha fc Zylstra (1991) 
suggested that these bursts are due to scattering of the X- 
ray emission by material evaporated from the accretion disc 
during the burst. This implies that double-peaked bursts can 
only take place over a certain range of inclination angles of 
the binary system. This model, however, does not explain 
the evolution of the temperature or the radius of the emit- 
ting area during the bursts, neither does it explain the fact 
that only a small fraction of the bursts are multi-peaked. 
iFisker et al.l l|2004') proposed that the double-peaked bursts 
are due to an interaction between the shell where the helium 
flash takes place, and a waiting point in the rp-processes of 
the shells above. This model can explain the multi-peaked 
bursts naturally in terms of a waiting point in the thermonu- 
clear reaction during the burst, but it is difficult to reproduce 
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the large dip o bserved between two peaks in some double- 
pea ked bursts ( Bhattacharvva fc Strohmaverll2006al ) . 

iBhattacharvva fc Strohmave 3 (|2006al lbh proposed that 
the double-peaked shape of the profile of some type-I X- 
ray bursts may be due to the latitude at which the ignition 
of the burst takes place: For bursts that ignite at a high 
latitude, the spreading of the flame on the surface of the 
neutron star is opposed by the flow of material accreted 
from the disc as this material moves towards the pole. As 
it approaches the equator, the burning front stalls momen- 
tarily, causing a minimum in the light curve. Eventually 
the burning front overcomes the opposing effect of the ac- 
creting material, and extends across the whole surface of 
the neutron star, leading to a double-peaked burst. This 
model can reproduce both the light curve as well as the 
spectral evolution of double-peaked bursts. It also provides 
an explanation of the millise c ond os cillations detected by 
IBhattacharvva fc Strohmaverl (|2006bl ) during the first peak 
of a double burst in 4U 1636-53. However, polar ignition (a 
crucial factor of this model) should occ ur at the highest ac- 
cretion rates (|Cooper fc Naravan|[2007l ). which makes it dif- 
ficult for this model to explain that the inferred mass accre- 
tion r ate at the time of the double-peaked b urst with oscilla- 
tions (IBhattacharvva fc Strohmaveij|2006b ) is lower than in 



19861). and all double-peaked bursts llGallowav et all l2008l : 
Maurer fc Wattsll2008l : IWatts fc MaureJl2007l . Zhang et al.. 



single -peaked bursts without o scillations (Watts fc Maurerl 
l2007f ). Notice, however, that ' Maurer &: Wattd ('2008') 
suggested that polar ignition can occur in a small 
rang e of lower accretion rates as well. Further- 
more, Im aurer fc Watts! I 20081 ) simulated burst light-curves 
for different ignition latitudes using a phenomenological 
model for the time-dependent surface-temperature profile 
o n the neutron-star surfa ce, and adopting the calculations 
oflSpitkovskv et al.l (I2OO2I ) for the speed of the burinng front. 
iMaurer fc WattT i 2008h also simulated light curves of bursts 
igniting quasi-simultaneously in two different places on the 
neutron-star surface, with delays between ignitions of 0.5 s 
to 4 s. None of their simulations produced a multi-peaked 
profile. 

Our detection of a triple-peaked burst in 4U 1636- 
53 casts doubt on the idea that the place of ignition of 
the burst is the cause of the mult i-peaked light curves 
(IBhattacharvva fc Strohmaverl [2006al ) . If this scenario ap- 
plies to the triple-peaked burst, the burning front must 
stall twice before it engulfs the whole surface of the neu- 
tron star. In their model of the double-peaked bursts, 
IBhattacharvva fc Strohmaverl (|2006al ) argued that the burn- 
ing front stalls as it approaches the neutron-star equator, 
and it regains speed after it crosses the equator. Clearly, a 
more complicated pattern of stalling is required in the case 
of a triple-peaked burst. 

Both double- and triple-peaked bursts in 4U 1636- 
53 have very similar properties: (i) The fluence of the 
triple-peaked burst reported here is 2.9 x lO"'^ erg 



the 



range 



X 10 erg cm (see 



Figure 4 in 

I Watts fc Maurer 2007). (ii) Except for the dou ble-peaked 
burst in IBhattacharvva fc Strohmaverl (|2006bl ). which is 
more or less a factor two brighter than the rest, all double- 
and triple-peaked bursts have similar peak-flux values, (iii) 
Neither the double- or triple-peaked burst are radius expan- 
sion bursts, (iv) Both triple-peaked burst s (the one reported 
here and the previous one reported in Ivan Paradiis et aU 



in prep.) occur when the source is more or less in the same 
region in the colour-colour diagram, accreting at relatively 
high inferred mass accretion rate. All these similarities sug- 
gest that the underlying mechanism for double- and triple- 
peaked bursts should be the same. Currently available mod- 
els are inadequate to explain the multiple-peaked bursts. 
It is clear that any model trying to explain the burst pro- 
files must also be able to explain the fact that only about a 
dozen out of ~ 250 bursts in the RXTE data of 4U 1636-53 
show multiple peaks, with only one out of 250 in the RXTE 
data, and in fact only two out of >4.5 x 10^ bursts observed 
since the beginning of X-ray astronomv (Ivan Paradiis et al.l 
1 19861 : 1 in't Zand et al. I I2OO4I : iGallowav et al.ll2008D . showing 
a triple-peaked profile. 
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